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Abstract. We report the high S/N observation on October 3, 2002 with XMM-Newton of the brightest X-ray 
flare detected so far from Sgr A* with a duration shorter than one hour (~ 2.7 ks). The light curve is almost 
symmetrical with respect to the peak flare, and no significant difference between the soft and hard X-ray range is 
detected. The overall flare spectrum is well represented by an absorbed power-law with a soft photon spectral index 
of r=2.5±0.3, and a peak 2-10 keV luminosity of 3.6tg j x lO^'^ ergs"\ i.e. a factor 160 higher than the Sgr A* 
quiescent value. No significant spectral change during the flare is observed. This X-ray flare is very different from 
other bright flares reported so far: it is much brighter and softer. The present accurate determination of the 
flare characteristics challenge the current interpretation of the physical processes occuring inside the very close 
environment of Sgr A* by bringing very strong constraints for the theoretical flare models. 
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' 1. Introduction 

■ Measurements of star motions in the central part of our 
, Galaxy demonstrated the presence of a high concentra- 
' tion of matter within a very compact region of 0.01 pc 
, at a distance of 8 kpc. These results prove the existence 
' of a central supermassive black hole of a mass of about 
, 3x10^ M0 at the dynamical center of our Galaxy (e.g., 
Schodel et al. 12002(1 . coincident with the compact radio 
source Sgr A*. Surprisingly this source is much fainter 
than expected from accretion onto a super-massive black 
hole. In particular in the 2-10 keV energy band its X- 
ray luminosity is only about 2.2x10^^ erg s~^ within a ra- 
dius of 1.5" (Baganoff et al. I2003al hereafter B03). This 
value may in fact be considered as an upper limit since 
this region contains other components such as stars (e.g., 
Eckart & Genzel lTM?! . hot gas, etc. Thus, Sgr A* radi- 
ates in X-rays at about 11 orders of magnitude less than its 
corresponding Eddington luminosity. Its bolometric lumi- 
nosity is only about 3x10"*^ LEdd (Meha & Falcke .20011 
Zhao et al. I2003|l . This has motivated the development 
of various radiatively inefficient accretion models to ex- 
plain the dimness of the Galactic Center black hole, e.g. 
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Advection-Dominated Accretion Flows (e.g., Narayan et 
al. [TM5|l . jet-disk models (e.g., Falcke & Markoff IMjUjl . 
Bondi-Hoyle with inner Keplerian flows (e.g., Melia et 
al. 12000(1 . The recent discovery of X-ray flares from Sgr 
A* has provided new exciting perspectives for the un- 
derstanding of the processes at work in the galactic nu- 
cleus. The first detection of such events was found with 
Chandra in October 2000. The flare has a duration of 
about lOks, with L(2-10keV)= l.OiO.l x lO^^ergs^i 
for the flare peak, i.e. about 45 times the quiescent state. 
(Baganoff et al. HUUTl hereafter BOl). The photon flare 
spectral index, obtained from the simultaneous fitting of 
the quiescent phases and the flare phase, was F=1.3^o'g 
(at 90% confidence level). A second significant flare was 
detected with XMM-Newton by Goldwurm et al. 1(2003(1 
who found a monotonic flux rise up to a factor of about 
20-30 in the last 900 s of the observation. They found also 
a rather hard spectrum with r=0.9±0.5 (at 68.5% confi- 
dence level), with a maximum observed 2-10 keV luminos- 
ity of about exlO^^'ergs-^ Baganoff et al. l|2003b(l esti- 
mated, with the analysis of several Chandra observations, 
that X-ray flares happen at a rate of 1.2±0.6 per day. 

We report here the detection and the analysis of the 
up-to-now brightest X-ray flare from Sgr A*, observed in 
October 2002 with XMM-Newton. The present data pro- 
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vide the highest signal to noise measurements of a Sgr A* 
X-ray flare, reported so far. 

2. XMM-Newton observation and results 

The Galactic Center was observed on October 3, 2002 
with Sgr A* at the center of the XMM- JVewton/EPIC field 
of view (see Turner et al. lMTTl and Striider et al. lMTTl for 
the description of the EPIC cameras MOS and FN, re- 
spectively). The MOS (exposure time ^16.8 ks) and PN 
(~13.9ks) cameras were both operated in the prime full 
window mode, with the medium and the thick filter se- 
lected, respectively. The data were processed with the 
XMM Science Analysis Software (version 5.4.1). X- 
ray events with pattern 0-12 and 0-4 are used for the 
MOS and PN, respectively. 

We find a strong increase of the count rate dur- 
ing about 3ks. To identify the counterpart of this X- 
ray flare, we improve the astrometry by cross-correlating 
the 8 PN sources detected in the 0.5-1.5 keV energy 
band during the whole observation with the 2MASS 
All-Sky Point Source Catalog using the eposcorr SAS 
task. We find a reliable 2MASS counterpart for each 
of these serendipitous PN sources, from which three 
are Tycho sources, by introducing boresight correction 
and rotational offset of (-1.2", -0.5"), and 2.7', respec- 
tively. The refined position of the X-ray flaring source is 
then aj2ooo= 17h45'"39.98% ^j2ooo= -29°00'28.4", with a 
90% confidence level error box of 2.3". The radio coun- 
terpart of Sgr A* at aj2ooo=17^45"40.0383'^ ± 0.0007^ 
'5j2ooo=-29°00'28.069"±0.014" (Yusef-Zadeh et al. 1999). 
is located only 0.8" away from the flaring source. Since 
no other X-ray source was found by Chandra in our error 
box, we conclude that the X-ray flare is very likely to be 
associated with Sgr A*. For the light curves and spectral 
analysis, we use an extraction region of 10" radius around 
Sgr A*, where the solar flare background is negligible with 
a maximum count rate of 6x 10""^ cts s~^. 



2.1. The EPIC light curve 

The EPIC light curve in the 2-10 keV energy range shows 
for the first 10 ks a quiescent period and then a very sharp 
and intense flare with a duration shorter than one hour 
(~2.7ks, see Fig. ^ top panel). The flare shape is almost 
symmetric relative to the peak. Both the rise and decay 
phases can be fitted by an exponential with characteristic 
times of 610 s and 770 s, respectively. We also observe a 
"micro-structure" at about 11,200 s with a significance of 
at least 3a on a time-scale as short as 200 s, which can 
further constrain the emitting region to a size of about 
7 Schwarzschild radii for the black hole mass of about 
SxIO'^Mq. The 2-10 keV count rate at the flare peak is 
1.74±0.13 ctss"^, i.e. 9.2^2 ^ times higher than the mean 
value measured during the quiescent period (0.19±0.04 
ctss^^). B03 using Chandra found that the inner 10" ra- 
dius region around Sgr A* contains six other X-ray point 
sources as well as a diffuse component, and estimated that 
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Fig. 1. EPIC light curves (MOS 1-l-MOS 2-hPN) and 
hardness ratio (HR) within a radius of 10" around Sgr A* 
position. The time binning is 100s, and the error bars in- 
dicate 1(7 uncertainties. Upper panel: The 2-10 keV light 
curve shows the quiescent and the flare periods. Middle 
panel: The 2-5 keV and 5-10 keV hght curves during the 
flare period. Lower panel: quiescent subtracted HR during 
the flare period. Long dashed line: mean HR value, and 
short dashed lines: la error bars. 



the flux from Sgr A* during its quiescent state represents 
about 10% of the total 2-10 keV flux. Therefore we in- 
fer for the present flare a count rate increase of a fac- 
tor of about 90 between the quiescent state of Sgr A* 
and the flare peak. The shapes of the light curves in 
the soft (2-5 keV) and hard (5-10 keV) energy bands are 
very similar (Fig. ^ middle panel). We report in Fig. ^ 
{lower panel), the hardness ratio (HR) during the flare pe- 
riod: HR= (hard— soft)/ (hard+soft). We removed the qui- 
escent level contribution, assuming that it is constant dur- 
ing the flare: mean count rates of 0.12±0.03 ctss^^, and 
of 0.07±0.03 ctss^^, respectively for the soft and hard 
energy bands. This corresponds to a HR value for the 
quiescent period of — 0.29±0.02. We find a flare HR of 
— O.liO.l, indicating a rather soft spectrum. No signifi- 
cant spectral evolution is observed between the [2-5] keV 
[5-10] keV ranges (a similar result is found for a different 
energy range, e.g. [2-3.5] keV). We find an indication that 
the flare presents a higher HR value, compared to the qui- 
escent, in a restricted time period of the rise phase (i.e 
11,000 < time < 11,400 s). However, this does not imply a 
change in spectrum of the flaring source since the emission 
within 10" during the quiescent period is not dominated 
by Sgr A*. 
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Table 1. Fit of the EPIC flare spectra, in the 2- 
lOkeV energy range, with absorbed power-law (pow), 
bremsstrahlung (brems), black-body (bb) or mekal mod- 
els, taking into account dust scattering. Frozen parameters 
are given between brackets, (a): Nh is units of 10^^ cm~^. 

(b) : Z is the metal abundance relative to the ISM one. 

(c) : The fluxes are unabsorbed and expressed in units of 
lO^^^erg cm~^ s~^. 
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xVd.o.f. 


pow 


2.0±0.3 
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104.3/123 


bb 


1.0±0.2 
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Fig. 2. EPIC spectra of the flare, within a radius of 10" 
around Sgr A* position. The spectrum for the quiescent 
period is used as background. A power-law model has been 
fitted taking into account ISM absorption and dust scat- 
tering. See Tabled for parameter fit values. 



2.2. Spectral analysis of the X-ray flare 

Fig. 121 shows the extracted spectra for the overall flare 
period. The overall flare spectra and the rise and decay 
phase spectra were binned to a minimum of 20 counts 
and 10 count per bin, respectively. XSPEC vl 1.2.0 is used 
for the spectral analysis. We use the response matrices 
from the rmfgen and ARFGEN SAS tasks. The errors 
and upper limits quoted correspond to 90% confidence 
level for one interesting parameter. The quiescent spec- 
trum including the contribution of Sgr A* in its quiescent 
phase and the contribution of the other sources in a radius 
of 10" around Sgr A* is subtracted from the flare spec- 
trum. Wc fit the data with an absorbed power-law model 
taking into account the dust scattering effects along the 
linc-of-sight, using the scatter model soon available in 
XSPEC (Predehl & Schmitt I1995|l . We fix for the scat- 
ter model, a visual extinction value Ay of 30 mag, as 
determined from IR observations of stars close to Sgr A* 
(e.g., Rieke et al. I1989;i . The interstellar medium (ISM) 
column density is fitted using the model tbabs/xspec 
(Wilms et al. 12000(1 . which includes updated cross-sections 
for X-ray absorption by the ISM. We also use their revised 
ISM abundances, which produce changes in ISM cross- 
sections of up to 30% with respect to solar abundances 
generally used as reference abundances for the ISM. This 
means that the absorption column values found here are 
about 30% higher than the ones obtained assuming so- 
lar abundances. The power-law continuum model with a 
soft photon index gives a good representation of the data 
(see Fig. 121 and Table and we find at the flare peak 
a 2-lOkeV luminosity of 3.6^q4x10'^^ ergs"^, assuming 
d=8kpc. Fixing r=1.3, as found in BOl, we obtain a non 
satisfactory data fit (Table . The presence of narrow 
FcKq, Gaussian emission lines ((T=10eV) are not statis- 
tically required (Ax^ <1 for one additional parameter). 
We find at 6.4 keV, 6.6 keV, 6.7 keV, and 7.0 keV upper 



limits on the equivalent widths of 33 eV, 122 eV, 123 eV, 
and 65 eV, respectively. Likewise the presence of a broad 
emission line near 6.7 keV ((T=250eV) and an absorption 
edge at 7.5 keV are not statistically required and we ob- 
tain EW<173eV and r <0.2, respectively. The data are 
also well fitted by a bremsstrahlung continuum model with 
a temperature of about 6.1 keV (Table 0. We also test 
an optically thin plasma using the mekal model with 
ISM abundances, letting the metal abundance vary as a 
free parameter, we find a good fit for a relative metal 
underabundance of less than 0.17. While fixing the rel- 
ative metal abundance to be equal or two times higher 
than the ISM abundances, the data fit is not satisfac- 
tory, with very high temperatures and therefore a large 
over-prediction of a H-likc iron line at 7 keV (Table ^ . 
Surprisingly, the spectrum can also be well fitted by a 
black-body with kT=1.47keV with no significant change 
of temperature all over the flare. The apparent size of 
the emitting area is about 10 km^ at the peak for a dis- 
tance of 8kpc. The spectral fits of the rise and decay 
phases of this flare are statistically consistent with each 
other with no change of the photon index {Mh fixed 
to 2x10^3 cm-2), with r=2.4±0.2 (xVd.o.f.=103. 0/122), 
and r=2.5±0.2 (xVd-o.f.=131.9/126), respectively. This 
confirms the lack of significant hardening according to the 
HR found. 

3. Discussion 

Several models have been proposed to explain the first 
Sgr A* X-ray flare observed with Chandra. For instance, 
it could be produced by a sudden enhancement of accre- 
tion through the inner hot circularized flow predicted by 
the model of Liu & Mclia ( 2002 ) . When the viscosity in- 
creases, synchrotron self-Compton (SSC) of the millimeter 
to sub-millimeter photons is the dominant process, and a 
soft spectrum is expected. Such flares would be coupled 
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with a similar increase in the mm and sub- mm ranges, 
which have not been yet observed (e.g., from March 2001 
to May 2002, Zhao et al. '20031. If the viscosity decreases 
then thermal Bremsstrahlung is the dominant process, and 
a hard spectrum would be generated, consistent with the 
slopes observed in the two other X-ray (much lower S/N) 
flares. Three possible flare mechanisms have been pro- 
posed by Markoff et al. H2001|l on the basis of a jet model: 
increased jet power or accretion rate, increased heating of 
quasi-thermal relativistic particles, or sudden shock accel- 
eration. The first scenario predicts a steep spectral index, 
and a huge radio flare for an accretion rate increase, and 
simultaneous flaring at all frequencies for an jet power in- 
crease, but this have never been yet observed. The second 
and third scenarii predict harder spectrum and were com- 
patible with the spectral index of the first Chandra X-ray 
flare. 

The flare luminosity excludes normal stars as the ori- 
gin of this event, but compact objects, although unlikely, 
cannot be totally excluded (see discussion in BOl). We 
rule out the possibility that the symmetrical shape of the 
flare light curve is produced by micro-lensing of a star lo- 
cated behind Sgr A*. Indeed from Zinnecker H2003|l 's for- 
mulae, we flnd that hour duration micro-lensing events 
imply very high tangential velocities for the lensed star 
and unphysical distances from the black hole, less than 
one Schwarzschild radius. For instance a velocity as high 
as 10,000 km s~^ corresponds to a distance of about only 
10^^ cm. The duration of the trailing edge of this flare is 
very short. If the physical process is thermal, this could 
be explained by the very rapid cooling of a plasma not 
magnetically confined, as proposed by Merloni & Fabian 

The total column density inferred from the spec- 
tral fit is about 1.7-2.0xl0^'^cm^^ for the power-law, 
bremsstrahlung and mekal models. The visual extinction 
Ay ^30 mag, inferred from IR observations of stars very 
close to Sgr A* (Rieke et al. I1989|l , corresponds to an ab- 
sorption column density along the line-of-sight of about 
GxlO^'^cm^^ (Predehl & Schmitt il995|ll . This discrepancy 
between Ay and Mh (also reported by Predehl & Triimper 
I1994|l may indicate an additional, probably related to 
the source, gaseous absorption from nearly neutral gas. 
A spectral analysis of other X-ray sources in the XMM- 
Newton field of view will enable to check whether this 
absorption excess is peculiar to Sgr A* (Predehl et al., 
in preparation). Nayakshin et al. I)2()()3|l proposed that 
such X-ray flares are due to the passage of stars through 
the accretion disk around Sgr A*, and that extra absorp- 
tion may be related to the neutral accretion disk pho- 
tosphere through which the star passes. The excess ab- 
sorption could also be explained by an overabundance of 
heavy metal elements from a putative Supernova remnant 
or He and other heavy elements from He stars located 
in the neighbourhood of Sgr A* (e.g., Eckart et al. 1995). 
However, for the blackbody model no extra absorption is 
required. 

We find also a soft X-ray flare spectrum with 



r=2.5±0.3, much softer than for the two other bright 
flares reported up to now (with lower S/N): r=1.3lo;g 
(BOl), and r=0.9±0.5 (Goldwurm et al. 2003). The 
present spectral index is consistent with the ones found 
during the quiescent periods of Sgr A* by Chandra: 
T=2.5toj in September 1999 (B03), and r=1.8l:[!;^ in 
October 2000 (BOl). Such a soft X-ray spectrum is pre- 
dicted by the scenario proposed by Liu & Melia ( 2002) in 
case of an increased viscosity of the accretion flow, and 
by the one proposed by Markoff et al. H2()()l(l in case of 
an increase of accretion rate or jet power. However the 
predictions of these two latter models are not consistent 
with observational constraints found at other wavelengths 
(e.g, Zhao et al. IMIHl Hornstein et al. Although, 
the Nayakshin et al. I)2()()3|l star-disk interaction model 
predicts a soft X-ray spectral index (P ~2.3) and a mod- 
erate Fe Kq, line EW of about 60 eV the allowed temper- 
ature values (kT=l and 70keV) are very flne-tuned (S. 
Nayakshin, private communication). Emission from a hot 
optically thin plasma with kT~5.7keV and a low relative 
metal abundance (<0.17 compared with the ISM abun- 
dances) is also consistent with the observed spectrum. But 
the high metal depletion is difficult to explain in view of 
the ambient, diffuse and metal rich emission and the pres- 
ence of close supernova remnants such as Sgr A East. If 
the flare is related to the accretion disk, the low metallic- 
ity could indicate a less evolved material in the disk. 

In conclusion, the present observation strongly con- 
strains the flare modeling of Sgr A*. 
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